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Abstract: The bottom pivot is a vital support device in the miter gate but often subject to poor 
lubrication and wear failures. Calculating the hydrodynamic lubrication characteristics of the bot- 
tom pivot is a complex three-dimensional (3D) problem, and most of existing models adopt simpli- 
fied assumptions to reduce the calculation difficulty. To solve this issue, this work develops a 3D 
model to calculate the hydrodynamic lubrication characteristics of the miter gate bottom pivot. The 
finite difference method is used to solve the oil film thickness and pressure distribution based on 
the spherical coordinates Reynolds equation. The component forces in three directions are calcu- 
lated from the pressure distribution and compared with the theoretical values to generate the cal- 
culation difference. Then, the genetic algorithm (GA) is used to minimize the difference to determine 
the optimal initial parameters for the 3D model. The analysis results show that the calculation ac- 
curacy can be significantly improved by using the optimal initial model parameters. When our ini- 
tial pressure is 5.64MPa, the results meet the engineering accuracy requirements. 
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Nomenclature 

h oil films thickness 

H dynamic viscosity of lubricating medium 
P oil film pressure 

U relative movement speed of the friction pair 
t time 

R distance from the point of the origin 

0, o spherical coordinates 

T shear stress 

Ey, Ey, E, eccentricity of the spherical bearing 

O, O,, Q, relative rotational angular velocity 

To ultimate dynamic shear stress 

C clearance of friction pair 
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1. Introduction 


Currently, a large number of miter gates suffer from structural damages due to ex- 
cessive bottom pivot wear [1]. As an important supporting component for the miter gate, 
the bottom pivot works in a harsh underwater environment and is prone to wear failures. 
Daniel [2] found excessive wear on the miter gate bottom pivot in the Dutch navigation 
locks while Zhao et al. [3] suggested that it is crucial to measure the surface wear of the 
bottom pivots to prevent unexpected failures and prolong the service life of the pivots. 
Xiaoxue Wu[4] improved the data set, changed the noisy data set into a clean data set, and 
then applied it to the security vulnerability report (SBR), and the prediction effect is much 
better than before.Sheng Chunfu[5] Taking ultra high molecular weight polyethylene 
(UHMWPE) fiber as the research object, they discussed the effects of wear frequency and 
yarn tension on its wear life.Ning Fanggang[6] proposed a method of using functional 
graphene / polyurethane composite coating (FG / PU) to improve the wear resistance of 
yarn to high modulus polyethylene fiber (HMPE) yarn. 

Due to unique structural features, the mathematical model of the bottom pivot can 
be simplified as a spherical bearing model. The lubrication characteristics of the bottom 
pivot can be analyzed by the spherical bearing lubrication. Pylos et al. [7] predicted the 
possible lubrication states of the wrist implants in the spherical bearing surface and the 
minimum oil film thickness under different working conditions using the elastohydrody- 
namic (EHD) lubrication theory. Deng et al. [8] used the Reynolds equation in spherical 
coordinates to inspect the lubrication of the bearing ball. Agrawal [9] added micro texture 
on the surface of the hybrid spherical thrust bearing and solved the modified Reynolds 
equation by finite element method. Wang and Sharma [10] solved the steady-state Reyn- 
olds equation on the spherical bearing using the finite difference method and relaxation 
iteration method. Kumar et al. [11] used the deterministic theory to analyze the influence 
of the surface roughness on the squeezed film lubrication of the spherical bearing. Huang 
et al. [12] studied the lubrication of the spherical crown pit texture model under different 
surface texture distribution parameters using the hydrodynamic lubrication theory. It can 
be learned from existing publications that the difference method is often used for the iter- 
ative solutions for calculating the Reynolds equation, resulting in low accuracy and heavy 
calculation time. The evolutionary algorithm can address this issue by properly optimiz- 
ing the initial model parameters in the difference method. Ji Y [13] combined the radial 
basis function with the NSGA-II genetic algorithm to carry out multi-objective optimiza- 
tion design for the turbine impeller. The efficiency and accuracy of the model has been 
improved. Tiwari et al. [14-17] used the genetic algorithm for optimum design of rolling 
bearings and spherical roller bearings. However, investigating the hydrodynamic lubri- 
cation characteristics of the bottom pivot using the genetic algorithm optimized finite dif- 
ference model has not been found yet . It is worth combining them to improve the calcu- 
lation performance on the hydrodynamic lubrication characteristics of the bottom pivot. 

This work develops a genetic algorithm optimized finite difference model to investi- 
gate the hydrodynamic lubrication characteristics of the bottom pivot. Firstly, the force on 
the bottom pivot is calculated through theoretical analysis. Then, the initial model param- 
eters are used to solve the Reynolds equation in the spherical coordinate, so as to obtain 
the film thickness and pressure distribution of the bottom pivot. Then the calculated re- 
sults are compared with theoretical values, and the difference is used as the feedback in- 
formation to optimize the initial parameters using the genetic algorithm. Changes in oil 
film and pressure of the pivot after parameters optimization are studied under different 
operation conditions. In this paper, a method combining finite difference method and ge- 
netic algorithm is proposed to solve the Reynolds equation. In the process of solving, there 
is no need to deal with the Reynolds equation dimensionless, but only need to replace the 
actual parameters. The model used can be applied to the solution of Reynolds equation in 


different cases, and the difference is less than 1.0x10°. 
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The remainders of this paper are organized as follows: theoretical model of the pivot 
is established in section 2 and the finite difference model is introduced in section 3. Sec- 
tion 4 presents the optimization procedure of the genetic algorithm. The numerical calcu- 
lation results are analyzed in Section 5. The main conclusions are drawn in section 6. 


2. Mathematical model of a bottom pivot 


The general form of the Reynolds equation is:[18] 


3 — 
vv )=v- ais a 
12u 


ôt 
(1) 
where h is the oil films thickness; 4 is the dynamic viscosity of lubricating medium; P 


is oil film pressure; U is the relative movement speed of the friction pair; t is the time. 


Replace X,Y and Zin the rectangular coordinates with 0 and É inthe spherical 

coordinates, it yields 

x=RsinOsin @ 

y= Rsin cos 

z=Rcosé 

(2) 

where R is the distance from the point of the origin;@ is the angle between the point 
and z-axis.? is the angle between the projection of this point in the XY plane and the Z- 
axis. 

The bottom pivot is lubricated with grease, and its lubricating medium can be con- 
sidered as a Bingham fluid, so the shear stress T can be expressed as 

T=T, +0 x 
* "06 (3) 
where To is the ultimate dynamic shear stress. 

Due to high viscosity of the lubricant grease and small gap between the bottom pivots, 
the grease flow can be considered as a laminar flow. The classical Reynolds equation can 
be derived from the Navier-Stokes equations if ignoring the curvature of the fluid film, 
inertial force and slip boundary. The Reynolds equation of the incompressible fluid in the 
spherical coordinates can be obtained as follows [18]: 
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where, Uo and U, are the linear velocity of the fluidin @ and Y directions, and R 
is the radius of the base pivot. 


Using Ex, y and £: to represent the eccentricity of the spherical bearing and @:, 


®, and ©®- to represent the relative rotational angular velocity, we can get 


h=c(1—é, sin @cos@—é, sin Asin g— €, cos 0) 


(5) 
—Ra* sing+ Ræ” cos ø 
2 (6) 


U- -Ræ* cos ġ cos 0 — Ræ” sin ġ cos 0 + Ra’ sind 
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where c to represent the clearance of friction pair, 


78 
79 
80 
81 


82 
83 


84 
85 
86 


87 
88 


89 
90 
91 
92 
93 
94 


95 


96 
97 
98 
99 
100 
101 


102 


103 
104 


105 


106 


107 


108 


109 
110 


Fy 


Figure 1. Brief structure drawing of miter gate 
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Figure 2. Schematic diagram of oil film thickness distribution (vertical view & front view) 


The miter gate structure is shown in Figure 1 and the oil film distribution is depicted 
in Figure 2. The bearing capacity of the bottom pivot is calculated as followings: 
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Based on Figure 2, Eq. (3) is processed and simplified to obtain Eq. (9): 
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For the convenience of calculation, Eq. (9) can be reduced as 
2 2 
A lop Lepp ar (10) 


P op o ô$ 


The coefficients in Eq. (10) are described as 
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3. Numerical Analysis 124 

The finite difference method is used to solve the Reynolds equation, and the differ- 125 
ence relationship is described in Figure 3. The solution area is divided into 90 meshes in 126 
the direction of @ and 360 meshes in the direction of ® . The step length along the O 127 


direction is 7R/180 and the step length along the ? directionis 7R/180. 128 
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Figure 3. Difference relationship 130 
131 
According to the difference relationship in Figure 3, the partial derivative of the pres- 132 
sure at point (i, i) can be expressed as: 133 
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00),, 240 
Op _ Pi ja T Pija 
ô$), 2A% 
(12) 134 
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Substituting Eq. (12) into Eq. (11), the relation between p of each node and p of adja- 135 
cent nodes can be obtained as 136 
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After the equation is discretized, Jacobi iteration is used to solve the differential equa- 
tion. 
k+l _ k ~k 
Pij =(1-a) pi +a Pi; (15) 


k k+l ~k 
Where Pij is the current pressure, Pij is anew pressure, Pi, j is the pressure ob- 


tained by Eq. (13), and is the decimal between 0 and 1. 

The boundary was determined by the following well-known scheme: 

if Pij <9, p50 

Adopts the discriminant method of relative convergence in the Reynolds Equation, 
we derive 
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where, n is the number of grids divided by rows and m is the number of grids divided by 
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columns, and € is a small constant; and in this work g=1.0x10*. 

The pressure on the bottom pivot edge is set to zero. In the initial condition, once the 
pressure is determined, the array of the NumPy library in Python is applied to calculating 
the component force of each pressure point by multiplying the elements, and then, sum 
the component forces in each direction as the total component force in this direction. 

Let us assume that there is only one load in the Z direction and only one rotational 
angular velocity around the Z-axis. At the beginning of the calculation, the film thickness 
distribution matrix is calculated by Eq. (5), and the obtained film thickness is substituted 
for Eq. (11) to calculate the values of A, B, C, D, E. Substitute the computed values into 
Eqs. (13) and (14) to calculate the pressure distribution matrix; repeat the iterative process 
of Eq. (13), and then, judge whether the iteration converges on Eq. (15). After convergence, 
a stable pressure distribution matrix is obtained. During the iterative process, in order to 
remove the negative pressure, the calculated negative pressure is reset to zero to meet the 
Reynolds boundary conditions [15] after several iterations in this study. The calculation 
procedure is shown in Figure 4. 
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Input Lubricant viscosity 
Radial clearance 
Load and rotational angular velocity, 
Rotational angular velocity 


At the beginning of the cycle, initialize 
the values of eccentricity ex, ey,ez. 


$ 


The initial parameters are plugged 
into the formula to calculate the 
film thickness distribution and then 
the pressure distribution. 


Genetic algorithm is used to 
optimize the initial 
parameters 


Determine whether it is 
balanced with external loads 


Stop 


Figure 4. Schematic of the solution process 


4. Model Initial Parameters Optimization 


Since the component forces calculated by the finite difference method is affected by 
the initial model parameters, such as the initial pressure and initial eccentricity, the Ge- 
netic algorithm (GA) is employed to optimize the initial parameters. Due to complex cou- 
pling effect, the relationship between the initial pressure and the force component is non- 
linear. Thus, a proper objective function is a must for the GA optimization. In this study 
the objective function is selected as: 


f(x)=min(f, — f.) (17) 


In the above formula, f, represents the force in the vertical direction and Í. rep- 


resents the force calculated by the finite difference method.A random number generating 
function is used to create 12-bit binary random data to represent as the GA initial "genes". 
If num is the actual number converted from binary to decimal, it will be mapped into a 


specified range by 


X = numx (X nax ` X min) + X min (18) 


min 


When iterating the difference equation, in order to obtain a good "gene", an evalua- 
tion information (i.e., the fitness function) is needed, which is expressed as 


8) = (F(X) = FO max) + 0.0001 (19) 


After the evaluation information is determined according to Eq. (19), possible solu- 
tions will be selected. The solutions closer to the theoretical value are more likely to be 
retained, and the corresponding probability of selection is higher. Parameter settings of 
the GA and the pivot are shown in the following tables. 


Table 1. GA settings 


Parameter Value 


Population 200 
Probability of mutation 0.9 
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The crossover probability 0.1 
Evolution generations 40 


Table 2. Bottom pivot material and iterative parameters 


Parameter Value 
Base pivot mushroom head material 40 Cr 
Bottom pivot tile material QT600-3 
Quality of gate 850 t 
Lubricant viscosity 645 Pa:s 
Elastic modulus of bearing bush 169 GPa 
Bearing bush Poisson's ratio 0.286 
Tensile strength of bearing bush >600 MPa 
Initial eccentricity 0,0.1,0.1 
The iteration error 0.0001 
Maximum search times 5000 


5. Results and discussion 


5.1. The initial pressure 


According to the above solving process, the author wrote a Python program to opti- 
mize several working conditions. Taking a specific type of gate as an example, the initial 
parameter Settings of simulation are shown in Table 1 and Table 2. Three-dimensional, 
two-dimensional, spherical and rectangular images of oil film thickness and pressure dis- 
tribution are shown as follows: 


Film thickness’mm Pressure’ MPa 
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(b) 
Figure 5. Three dimensional schematic diagram of oil film thickness and pressure 
distribution 
Figure 5 shows that the oil film thickness gradually decreases to the @ direction. 
Since the initial parameters of the are preset, the oil film thickness along the P direction 
will first decrease and then increase. 


The normalized results of the maximum stress are shown in Figure 6. 
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Figure 6. Schematic diagram of the changing trend of normalized pressure and film thickness 

According to Eq. (13), during each iteration, the pressure at the center point in the 
grid plane is affected by the surrounding points. Therefore, different initial pressures will 
affect the calculation results. Other parameters are shown in Table 2. Tet the initial pres- 
sure from 1 MPa, 2 MPa, 3 MPa... to 10 MPa, and observe the impact of different initial 
pressures on the results. The results are shown in the following figure: 
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The initial pressure 


Figure 7. The effect of the initial pressure on the result of the force-splitting calculation 

As can be seen from Figure 7, the initial pressure and the relationship between com- 
ponents is a linear approximation, the increase of initial pressure to the element of differ- 
ent direction is different. Therefore, when calculating the bearing capacity of the bottom 
pivot, choosing an appropriate initial pressure is conducive to improve the calculation 
speed. 


5.2. Eccentricity 


In the actual process, the bottom pivot is affected by external load and its gravity 
resulting in different eccentricities. That the lubricating oil film produces a wedge effect, 
affecting the regular operation of the bottom pivot. According to Equation (4), the oil film 
thickness at any point in the bottom pivot is related to its eccentricity. To analyze the in- 
fluence of film thickness on the lubrication performance of the bottom pivot, different ec- 
centricity ratios are selected to observe the changes of force components in each direction, 
as shown in the figure below: 
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Figure 9. Influence of eccentricity in the Y direction on component force under different initial pres- 239 
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Eccentricity in the Z direction 


(c) 
Figure 10. Influence of eccentricity in the Z direction on component force under different initial 
pressures 


It can be seen from Figure 8, 9, 10 that when the eccentricity in the X direction in- 
creases, the bearing capacity in the X and Z directions increases and the bearing capacity 
in the Y direction decreases. The increase in eccentricity in the Y direction will also lead 
to the increase in bearing capacity and the greater the initial pressure, the more obvious 
the impact on bearing capacity. When the eccentricity in the Z direction increases, the 
bearing capacity of the X direction increases and the bearing capacity of Y and Z directions 
decreases. 


5.3. Rotational angular velocity 


When the spherical bearing does not rotate, there is only a static pressure effect. 
When it turns, the dynamic pressure effect is produced due to the existence of a wedge- 
shaped oil film, which affects the bearing capacity. Under the condition that other param- 
eters remain unchanged, the rotating speeds are 3° / s, 6° /s,9° /s, 12° /s and 15°/s. 
the bearing capacity is shown in the figure below. 
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(c) 
Figure 11. Influence of rotational speed on component force under different initial pressure 
As shown from Figure 11, with the increase of rotating speed, the bearing capacity in 
the X direction also increases, and the greater the initial pressure, the more obvious the 
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improvement. But the Y direction is the opposite. When the initial pressure in the Z direc- 
tion is less than 4MPa, the higher the speed is, the greater the bearing capacity is, but it 
will decrease when the initial pressure is greater than 4MPa. 


5.4, Optimization process analysis 


Even if the pressure obtained by the difference method converges, the obtained so- 
lution is not necessarily the fundamental solution to the problem, so the feedback pro- 
cesses needs to be used. The parameters of the genetic algorithm are 200 population sizes, 
40 number of evolutions, and 0.6 and 0.9 cross probabilities, respectively, and the goal is 
to minimize the difference between the carrying capacity and the theoretical value of the 
Z direction. The variation process of the initial pressure in the difference method of the 
evolution times is shown in the following figure: 


—— The crossover probability is 0.6 
The crossover probability is 0.9 
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Initial pressure (MPa) 


Evolution times 
Figure 13. Relationship between evolution times and initial pressure 
It can be seen from Figure 13. The initial pressure obtained when the crossover prob- 
ability is 0.9 gradually tends to be stable after 30 evolutions. In contrast when the crosso- 
ver probability is 0.6, it will not be stable after 30 changes. Therefore, this paper adopts 
the crossover probability of 0.9 for optimization. 
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Figure 14. Comparison diagram of optimization results 

As can be seen from the above figure, the difference between the calculated value 
and the theoretical value is quite different to different schemes. Schemes 1 and 2 are the 
results obtained after genetic algorithm optimization, and schemes 3 to 10 are the results 
of manually setting the initial value. It can be seen from the figure that the genetic algo- 
rithm has higher accuracy, and there is no need to try it many times manually. After GA 
optimization, when the initial pressure is about 5.64MPa, the ratio of the difference be- 
tween the estimated value and the hypothetical value of the hypothetical value is less than 


1.0x10°. 


6. Conclusions 


A numerical calculation model of hydrodynamic lubrication characteristics of miter 
gate bottom pivotis proposed to this paper. The effects of initial pressure, eccentricity and 
rotating speed on oil film thickness and pressure distribution are analyzed. In the iterative 
process, the increase in initial pressure will significantly improve the bearing capacity for 
the z-axis direction, and the increase in eccentricity in the X and Y directions will increase 
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the component forces along these directions. The greater the initial pressure, the faster the 
component increases. At the same time, the increase in Z-direction eccentricity will reduce 
the bearing capacity. In addition, the component forces in the X and Z directions will in- 
crease as the rotation increases. But the elements in the Y direction will decrease. Finally, 
the genetic algorithm is used to optimize the initial pressure, so that the error between the 
theoretical value and the calculated value is significantly reduced, and the calculated re- 
sults and external load meet the engineering accuracy requirements. Through the calcu- 
lation of finite difference method and genetic algorithm, the oil film thickness distribution 
and pressure distribution are consistent with the expectation. Because the load on the bot- 
tom pivot is large and vertical, the film thickness at the center of the ball is the thinnest 
and the pressure in this area is the largest.We can use more effective optimization algo- 
rithms. When there are more data, we can use machine learning and deep learning to train 
larger models to improve the accuracy and adaptability of the models. 
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